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Abstract 



We present a detailed analysis for the Langevin dynamics of a spherical spin- 
glass model (the spherical Sherrington- Kirkpatrick model). All the spins in 
the system are coupled by pairs via a random interaction matrix taken from 
the Gaussian ensemble. 

One finds that for a general initial configuration the system never reaches 
an equilibrium state and the theorems associated to 'equilibrium dynamics' 
are violated. Only very particular initial conditions drive the system to equi- 
librium. 

The weak ergodicity breaking scenario is demonstrated for general 'non- 
equilibrium' initial conditions. Two-time quantities such as the auto- 
correlation function explicitly depend on both times. When the time dif- 
ference is short compared to the smaller time one finds stationary dynamics 
with time translation invariance and the fiuctuation-dissipation theorem sat- 
isfied. Instead, when the time difference is of the order of the smaller time one 
finds non-stationary dynamics with aging phenomena, the system remembers 
the time spent after the initial time (the quench below the critical tempera- 
ture). Interestingly enough the short time- difference dynamics (FDT regime) 
for non-equilibrium initial conditions is identical to the relaxation within the 
equilibrium states obtained for the particular 'equilibrium' initial conditions. 
This points to a self-similaririty of the energy landascape. 

In addition we analyse the effect of temperature variations on the be- 
haviour of the correlation function and energy density. In this way we are 
able to make some comparisons between experimentally observed results and 
exact calculations for the model. Somewhat surprisingly, this simple model 
captures some of the effects seen in laboratory spin-glasses. 
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I. INTRODUCTION 



Spin-glasses, like many other complex systems^, are essentially out of equilibrium on 
experimental scales. Below the critical temperature experiments show 'aging effects', or the 
dependence of the response of the samples on their history since the temperature quench^'"^. 
The simplest way to observe aging phenomena is through the zero-field cooling (ZFC) ex- 
periment in which the sample is cooled in zero field to a sub-critical temperature at time to'-, 
after a waiting time a small constant magnetic field is applied and subsequently the time- 
dependence of the magnetisation is recorded. Its 'mirror' counterpart consists of cooling the 
sample at the the quenching time to = in the presence of a small constant magnetic field, 
keeping the sample in the field up to the waiting time and then measuring the decay 
of the thermoremanent magnetisation (TRM). In both cases, the 'older' the systems (the 
longer the waiting times) the slower their relaxations: the systems age^'"^. 

Other, more sophisticated, experiments include variations of the temperature, mainly 
temperature jumps and temperature cyclings below the critical temperature Tc during the 
waiting time. In the temperature jump experiments'*'^''' the samples are kept at a constant 
temperature T during a waiting interval [0, tj]. At tj the temperature is changed to T -|- ST. 
The sample is still kept waiting at this constant temperature during the time interval [tj, t^]. 
Finally, at the magnetic field is switched on or off and the ZFC or TRM decays are 
measured, respectively, at the final temperature T -\- ST. In the simplest setting'* tj = t^. 

In the temperature cycling experiments, a temperature cycle is performed during the 
total waiting time; T[t) is T for times between to and t^,!, T -\- ST for times between t^i 
and and again T for all subsequent times. The magnetic field is switched at > 
according to the ZFC or TRM settings^'^. 

In all these experiments, the notion is to find an effective waiting time corresponding 
to a system that has undergone temperature variations. This effective waiting time can be 
estimated by comparing the decays of quantities such as TRM with the same system kept 
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at constant measuring temperature. Whether to find the efi'ective waiting time is possible 
depends on the precise protocol of the experiment"^. 

The ZFC and TRM experiments'*'^ with temperature jumps at tj = unequivocably 
show that waiting at a higher temperature [ST < 0) disfavours aging and makes the materials 
respond as younger ones while waiting at a lower temperature [ST > 0) favours aging 
and makes the materials respond as older ones'* as to compared to a system kept at a 
constant measurement temperature T + ST. The response depends on the magnitude and 
duration of the heat pulses. The position of the maximum of the logarithmic derivative of 
the magnetisation decay is sometimes associated to an effective age of the system. After 
temperature jumps one observes a displacement of this maximum towards smaller (bigger) 
times^ when ST > [ST < 0). 

Grosso modo this is also observed by comparing the full TRM decays with and without 
jumps^'''. However, a more detailed analysis of the TRM decay curves after temperature 
variations shows more subtle effects and provokes controversy since there is no agreement on 
its dependence on the sign of the temperature variation. The Uppsala group'*'^ claims that 
their measurements of the rate of change of the decay of the magnetisation are symmetric 
with respect to the sign of ST. On the other hand, the Saclay-UCLA group studies the 
decay of the magnetisation and claims that it has an explicit dependence on the sign of ST^ 
i.e. that the response is asymmetric^'''. In Refs. 6 and 7 it is concluded that a short, as 
compared to the initial waiting time, heat pulse [ST > 0) partialy reinitializes aging while a 
negative temperature cycling [ST < 0) freezes the system into the state reached during the 
initial aging process. In other words, to lower the temperature is said to be equivalent to 
a, maybe partial, new quench. Instead, raising the temperature does not produce such an 
effect, (see Ref. 3 for a more detailed discussion on this discrepancy). 

The out of equilibrium dynamics of the standard finite dimensional model for spin-glasses, 
viz. the 3-dimensional Edwards- Anderson model (3D-EA), has been extensively studied 
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numerically® Simulations for the TRM decay and for the two-time auto-correlation func- 
tion decay at constant temperature®'^ and after temperature changes during the waiting 
time^°'^^ have been performed. In addition, simulations of the dynamics of the hypercu- 
bic spin-glass cell - a mean-field model for large dimensionalities - with a very schematic 
discussion of the effect of temperature shifts at t^, have also been carried out^^. 

Numerically, there is again agreement on the infiuence of the magnitude and duration 
of the heat pulses on the response of the 3D-EA, but there are different opinions as regards 
the symmetry or asymmetry in the response. In Ref. 10 it is claimed that the numerical 
results for the rate of change of the magnetisation® are symmetric. Conversely, in Ref. II 
numerical support for asymmetries in the response is given. 

The aging effects and in particular the effects of temperature variations have been inter- 
preted with various phenomenological models. 

The Uppsala group interpreted their experimental results along the lines of droplet 
models^"^^^^, i.e. in terms of time-dependent domain growth with the added assumption of 
chaoticity. These models predict a symmetric response to temperature changes. The Saclay- 
UCLA group interpreted their asymmetric observations with a 'hierarchical' modeP inspired 
by the replica solution of mean-field spin-glass models^^. In this model it is assumed that 
the system relaxes in a rough energy landscape with an hierarchical organisation of many 
metastable states. 

A semi-phenomenological approach to spin-glass dynamics based on the Parisi solution 
for the Sherrington-Kirkpatrick model has been proposed^''. One constructs a Markov chain 
on a hierarchical tree with jump rates which are quenched random variables chosen from 
a Levy distribution; the index of the distribution depending on the level to which the 
jump takes place. One associates an overlap between states depending on the ultrametric 
distance between them and then computes the correlation function within this framework. 
This model clearly exhibits aging phenomena and a simple two-level tree (corresponding to 
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two step replica symmetry breaking) may be used to account for many of the experimental 
observations. Temperature changes are implemented by a reorganisation of the tree structure 
via changes in the Levy law indices and overlap variables. 

Only recently attention has been paid to the analysis of the off-equilibrium dynamics of 
mean-field spin-glass models. The long-time analytic solution of the dynamic equations of 
the p-spin spherical modeP®, for p > 3, and for the Sherrington-Kirkpatrick modeP^, have 
been worked-out. In addition a detailed numerical analysis of the dynamical equations for 
a particle moving in an infinite dimensional random potential with long-range correlations 
has been carried out^° (see also Ref. 21). In these studies the dynamical equations are 
those associated to the relaxation of the system at constant temperature starting from a 
given initial condition. The effects of temperature variations have not yet been studied with 
mean-field spin-glass models. 

In this paper we shall study the dynamics of the spherical Sherrington-Kirkpatrick model 
or the p = 2 spherical spin-glass modeP^'^"^ for general initial conditions. The model is simple 
enough to be solved exactly statically and dynamically and allows us to examine in a simple 
analytical way some of the experimental scenarios. 

From the static point of view the model is extremely simple. Its energy has only two 
minima, corresponding to the configurations with maximum (minimum) projection on the 
direction of the eigenvector associated to the maximum eigenvalue of the (Gaussian) in- 
teraction matrix. Within the replica formalism the model is solved exactly with a replica 
symmetric ansatz^^; a result that would suggest its triviality from the dynamical point of 
view. Despite this, the model's dynamics is extremely interesting. For almost any initial 
condition it is out of equilibrium, and the evolution does not lead it to equilibrium. The 
model exhibits aging effects in the two-time auto-correlation function, i.e. even asymptoti- 
cally it depends explicitly on the waiting time. Depending on the initial conditions, rather 
than the exponential types of decay in the energy-density, correlation functions, etc. one 
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observes power-law decays. This is because the characteristic equihbriation time of the sys- 
tem (for such initial conditions) is infinite. Such an exponential decay with rate given by the 
inverse of this equilibriation time therefore cannot appear; the only other time-scale that 
can appear in correlation functions is the waiting time. 

The aim of this paper is two-fold. On the one hand we probe some of the assumptions 
used to obtain analytical results for the long-time dynamics of more complicated mean-field 
spin-glass models^®'^^; namely, the weak-ergodicity breaking hypothesis^''. On the other 
hand, we study the effects on this simple model of temperature jumps and temperature 
cyclings during the waiting time and we compare these analytical results with the associated 
experimental and numerical measurements. We do this via analysis of the decay of quantities 
such as the energy density and the auto-correlation function. 

This being the 'simplest' mean-field spin-glass model, it is not expected to reproduce the 
experimental behaviour in every detail. Surprisingly enough, some of the strange features 
of real spin-glass behaviour are captured by it. 

The paper is organized as follows. In Section II we present the model and describe 
the quantities of interest. In Section III we present the results for constant temperature 
T[t) = T . In Section IV we study the effects of temperature variations on the asymptotic 
decays. Finally, in Section V we present our conclusions. 



7 



II. THE MODEL 



The p = 2 spherical model is defined by the Hamiltonian 

i7 = J.jS.Sj , (II.l) 

where s^, i = 1, . . . , are the spherical spin variables constrained such that J2iLi -^1 = 
The couplings between the spins are given by the quenched random variables Jiy 

To study the dynamical evolution it is convenient to diagonalize the coupling matrix Ji^ 
and to work with the time-dependent projections of the spin configuration sit) = {siit)} 
onto the J-eigenvectors = // • sit). 

The dynamics for the model is defined via the usual Langevin equation, which when 
projected onto the eigenvectors // becomes 

ds,{t) 



dt 



{ji-z{t))s,{t) + h,{t) + Ut)- (11.2) 



/i is the eigenvalue associated to the eigenvector hf^it) represents an external magnetic 
field, zit) is a Lagrange multiplier enforcing the spherical constraint and (^^(t) is the thermal 
noise with zero mean and correlation given by 

{Ut)Ut')) = 'iT{t)8,,8{t-t'). (II.3) 

Tit) is the (possibly time-dependent) temperature. Hereafter we shall use ( • ) to represent 
the average over the thermal noise. 
The general solution to Eq. (II. 2) is 



(t) = s^{to) e^(*-*°) e-^'l + f dt" e^(*-*") e" '^"'"("'^ {h^{t") + i^{t")) . (II.4) 



The two-time auto-correlation function is defined as usual: 
C{t,t') = (1/iV) ^Y^iLi{^i{t)si{t'))^ ^, where the square bracket indicates averaging over 
the disorder Jij. It can be expressed in terms of the eigenvalues of J as: C{t^t') = 
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/ dfi p[fi) (s^(t)s^(t')), with p[fi) the eigenvalue density. In the absence of an external 
field, using Eq.(II.4) one obtains 



(( is^ito)Y e^(*+*'-2*°) )) + 2 r"^*'* ^ dt" Tit") Tit") (( e^(*+*'-2*") )) 



(11.5) 



where (( • )) stands for / dppifi) ■ and 

Tit) = exp (2 ^* dt" zit")^ (II.6) 

may be computed self-consistently from the spherical constraint C(t,t) = 1, as the solution 
to the following Volterra equation of the second type: 



Tit) = (( i-s^ito))^ e2^(*-*°) )) + 2 / dt" Tit") Tit") (( e^^^*"*") )) . (II.7) 



This immediately implies Tito) = 1. In what follows we shall take the initial (quench) time 
to be zero: to = 0. 

In this paper we shall analyse the case in which J is a symmetric matrix with elements 
which are independently distributed Gaussian random variables with zero mean and variance 
proportional to this choice gives the model a well defined thermodynamic limit. The 

probability distribution function for the eigenvalues is then given by the Wigner semi-circle 
law24 



K/i) = :^^^ /iG[-2,2]. (II.8) 

However, as far as possible, we shall keep a general distribution pip). This will ultimately 
be useful in studying the effects of different, e.g. wide^^ distributions^^. 

We study the long-time dynamics by calculating the long-time behaviour of 
• The time dependence of the energy density and of its relaxation rate. In the absence 
of an external magnetic field one can show that the energy density is given by 
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m = \ im - z{t)) = \ [rii) - i I in(r(t))] (11.9) 



and its relaxation rate is given by 



|M = ! f^r«)-i|li„(r«))V (luo) 

dint 2 \dt ^ ' 2 dt^ ^ ^ ^ ' 



The (averaged over the noise) staggered magnetisation 



{s,{t)) = '-^^ . (11.11) 

/r(t) 



• The two-time correlation function (II. 5). We shall be interested in the aging behaviour 
contained in C(t,t') for t and t' large. Choosing t>t\ a, useful expression for C is the 
following: 

1 



C{t,t') 



r(i±^)-2/ ' Jt"r(t")r(t") ((e^(*+*'-'*"))) 



11.12) 



We shall also study the decay from the initial conditions characterised by C(t,0). 

• The equal-time overlap between two real replicas and cr^ that evolve with the same 
thermal noise up to a 'waiting time' t^,, after which they evolve in two independent 
realisations of the thermal noise, say S^it -\- 1^) and ^[t -\- 1^): 

Q{t + t^,t + t^) = {{ {s^{t + t^) a^{t + t^))^-^ )) 

1 

(( s^{o) a^{o) e2^(*+*») + 2 j^^ dt" T{t") ^Jv s{t")T,{t") g^Mw^-*"))) 

(11.13) 



The realisation of disorder is the same, i.e. the interaction matrix Jij is the same for 
both replicas. 
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• The response function at time t to the perturbation by a small magnetic field applied 
at time t' 



Sh.it'} 



h=0 



(11.14) 



h=0 



hf_i{t) = fi ■ h[t). Taking the variation of {s^[t)) with respect to h^[t') from Eq.(II.4) 
one obtains 



G{t,t') 



1 



\ r.(t) 



2 (r„(t))3/2 6h,{t') 



3,(0)e-* )) . (11.15) 



h=0 



The presence of a field modifies the function F. However, explicitly computing the 
functional derivative 6Tfi[t)/6h^[t') from the defining equation and taking the ther- 
modynamic limit, one gets ST fi[t) / 6h^[t')\fi=o = 0. Hence, 



((e^(*-*'))) 



(n.i6) 



the subindex o indicating the absence of an external field. 



The decay of the TRM (assuming a suitably well behaved response function) is then 
given by 



Jo 



(11.17) 



All these functions contain information about the nature of the dynamics and about the 
'geometry' of the energy landscape of the model. 
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III. CONSTANT TEMPERATURE 



In this section we shall analyse the evolution of the system at constant temperature. For 
the initial conditions Si(0) at the quenching time to = we choose two configurations that 
are representative of the different possible behaviours of the system: 

• 'Uniform' initial condition. 

s^(0) = 1, V/i =^ Non-equilibrium dynamics (IH-l) 

• 'Staggered' initial condition. 



g^d — i^^j I Non-equilibrium dynamics for a ^ 2 

Equilibrium dynamics for a = 2 



The 'uniform' initial condition that has a constant an equal to one projection onto each 
eigenvector of the interaction matrix J is a random initial condition when written in the 
original basis. It corresponds then to the 'realistic experimental' initial condition from which 
the samples evolve after the rapid quench at to- 

The strange form of the initial conditions in the staggered case follow from the re- 
quirement that the initial conditions satisfy the spherical constraint. In order to have an 
equilibrium dynamics, the initial condition must have a macroscopic condensation onto the 
maximum eigenvalue (// = 2 for p[fi) given by Eq.(II.8)). We shall show that only this par- 
ticular initial condition leads the system to equilibrium, by studying the initial conditions 
needed to have a time-homogeneous relaxation of the correlation and response functions, 
two typical features of the equilibrium dynamics. 

The function T[t) can be obtained for general initial conditions using Laplace transform 
techniques to solve Eq. (II. 7). We derive these results in Appendix A. 
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A. Uniform initial condition 



The function r(t), for all times t and T < Tc = 1, is given by 



T(t) = -y k T'^ . (IIL2) 

k=0 



For large times the energy density is 

Sit) ~ i r - 1 + 1^ (III.3) 

and hence asymptotically it tends to the equilibrium value^^ Seq = | T — 1 with a power law 
decay. 

The staggered magnetisation given by Eq. (11.11) behaves, asymptotically, as 

{s,it)) ~ iiTr'f' qEA (2t)3/\(^-2)* , (IIL4) 

T < Tc and qEA = ~ T) is the Edwards- Anderson order parameter^^. If // 7^ 2, it decays 
exponentially with time. Instead, if // = 2, the staggered magnetisation associated to the 
maximum eigenvalue grows with time as a power law: {s2{t)) ~ t^^'^. The system condenses, 
on average, onto the maximum eigenvalue. 

The rate of decay from the initial conditions is characterised by 

C(t,0)~ Q'^' qEAt-'^\ (III.5) 

i.e. a power law decay of the correlations between the system and its initial conditions - 
this is typical of non-equilibrium dynamics. 



The correlation function (II. 5) is given by 

'h{2{t + t')) , rt' l,[2{t + t' - 2r)) 



and asymptotically (large time t') it behaves as 
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+ 2T [ drViT 
Jo 



t + t' Jo ^ t + t' -2t 



(111.6) 



C(t,t')~2y2^^^^ [l-rinti(2t(l-A),2t(l + A))] , (IIL7) 
where A = t' jt and 

It is clear Irom Eq.(III.7) that even at zero temperature the correlation function shows 
an 'aging' behaviour: C(t,t') = C(A) and hence C(t + t^^t^) = C{t/tw). It can be also 
proven that dC{t,t')/dt < and dC{t,t')/dt' > 0, for all times t,t' such that t > t'. These 
are two of the 'weak-ergocicity' breaking properties''''^®. 

For large t there are three t'-regimes: 

• Finite t' or A = t'/t 

C(t,0~/(r,0r3/4 (IIL9) 



wi 



ith 



/(r,t') = (47r)-'/^ qEA 2^/^ 



2t' 



1 + 2T r dt"T{t")e-^'" 
Jo 



[111.10] 



Having already taken the limit t — )• oo one can now consider the limit t' — )• oo; thus 
C{t,t') ~ 2V2 qEA ^^^y^, ~ 2V2qEA y?'^ (III.ll) 



The correlation function decays to zero with a power law. Hence, given any finite time 
t' there exists a big enough subsequent time t such that the correlation function has 
decayed to zero. 

• Large t', it — t')/t ^ 1, A = 1 — (^7^) — ^ 1; i-e. t and t' are relatively close to each 



other. 



Cit, f] ^C{t-t') = l-T duj e-" . (III.12) 

Jo (jJ 
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Up to zeroth order in [t — t')/t^ tlie correlation function satisfies time-translation 
invariance (TTI). This is the so-called FDT (fluctuation dissipation theorem)-regime^®. 
If t = t', A = 1 and C(t,t) = 1, as expected from the spherical constraint. If we take 
the limit (t — t') — )• oo - determining the end of the FDT-scale - we obtain 

C{t-t')r^qEA (III.13) 
as can be clearly seen in Fig. 1. 

This property completes the weak-ergodicity breaking scenario^'''^^. That is to say 

linir^oo C{t + tiu,tiu) = V tu, finite 

liint^^^C{T + t^,t^) = qEA V T finite. 

The initial rate of decay of the correlation function is 

li„,^ = -T. (III.14) 

t'^t dt ^ ' 

The usual fiuctuation-dissipation theorem is satisfied for this range of times as we shall 
show below. 

• Large t', it — t')/t finite and A finite, i.e. widely separated times t and t' . 

If A < 1 the full Fq.(III.6) holds. For large t, inti(2t(l - A), 2t{l + A)) is a function of 
A. Hence, for T < Tc = 1 the correlation function depends explicitly on t^. This is the 
'aging' regime. The rates of decay dC{t^t') / dt and dC{t^t') / dt' are proportional to 
— l/t and 1/t, respectively, and hence the auto-correlation has a very slow variation. 
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Fig. 1 C(t + tyj^tyj) VS. T in a log-log scale, tyj = 30, 100,300, 1000 (from left to right) at 
T = 0.3. qEA{T = 0.3) = .7. 

In Fig. 1 we present a plot of C(t + t^^t^) vs. r in a log-log scale, for constant 
temperature T = 0.3 and waiting times = 30,100,300,1000,3000. The curves are the 
typical aging-curves already observed in the Monte Carlo simulations of the 3D Fdwards- 
Anderson modeP and the D-dimensional hypercubic spin-glass celP^. r is the time elapsed 
after the 'waiting time' t^. The waiting times in these figures are enough to show the 
asymptotic behaviour, the same behaviour is reproduced for longer waiting times. In all 
the curves two clearly different time regimes appear: for times r much smaller than the 
auto-correlation has a fast decay from C = 1 to C ~ qEA = 1 — T = .7, while for r > t^, it 
has a slow decay from C ~ qEA to zero. The waiting time acts as the time-scale determining 
the length of the plateau in Fig. 1, i.e. the length of the FDT regime. 



To compute the overlap between two (real) replicas we make them start from the same 
initial condition, the uniform one, and let them evolve in the way described in Section II. 
Then F^ = T^j and 

1 



Q{T + t^,T + t^) = 1-2T 



^ ' 2(T + U-t") 



(111.15) 
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If tu, = then for large times r, Eq. (III. 15) implies: 

g(T,T)~9|^. (111.16) 

For a non-zero t^, the behaviour depends on t^,. If t^, is large one can show that there are 
two relevant r-regimes: 

Ifl<T<t^ Q(T + t^,T + t^) ~ C{T + t^,t^) 

If T > > 1 Q{t +t^,T + tyj) ~ > C(t + tyj.tyj) 

Fig. 2 shows, in a log-log scale, the Q decay for t^, = 0, 10, 30, 100 at constant temperature 
T = 0.6. One can there see how the r — )• oo limit of the curves depends on t^,, and how it 
approaches qEA for increasing t^,. 



Q 




0.01 0.1 1 10 100 

t 



Fig. 2. The overlap Q{t -\- t^^r -\- t^) vs. r in a log-log scale at constant temperature 
T = 0.6 for = 0, 10, 30, 100 (from bottom to top). qEA{T = 0.6) = 0.4. 

These results go in the same direction as those obtained numerically for the SK modeP'' 
and analytically for the 0{N) modeP^. The systems escape from themselves faster than 
they do from each other. However, in this model, the asymptotic value of the overlap Q 
at large t^,, when r — )• oo is qEA while for the SK modeP'', Q{t -\- t^^T -\- t^) ~ when 
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the same limits are considered. For the p = 2 spherical model this implies that there is a 
long-term drift in the energy landscape which causes the two replicas to follow roughly the 
same route. (One could imagine the replicas moving through channels or through the same 
sets of valleys in the energy landscape.) This demonstrates the relative simplicity of the 
energy landscape in the sense that a more chaotic landscape {e.g. that of the SK model) 
would possess a more complicated set of local-free energy minima which allow the replicas 
to become more widely separated asymptotically. 



The response function reads 

rui')- [m h{2{t-t')) 

For large times t > t', the r.h.s. of (IL16) is independent of the temperature and 



;iIL17) 



:iiLi8) 



f t-f 

One can demonstrate from this equation that: 

• If t' is finite the response function decays as the power law G[t^t') oc t~^^'^. 

• In the FDT scale (t - t')/t ^ 0, A ^ 1, 

Git, f) ~ e-2(*-*') ^'^^^^Sf^^ , (111.19) 

i.e. it satisfies time-translation invariance. Furthermore, from Fqs.(III.12) and (III. 19) 
one can immediately see that the FDT is also satisfied in this scale: 

T G{t - f) = ^^^^ . (III.20) 

• Finally, in the aging regime, 

G(M') ~ (Y^^ ^-^/^ (III.21) 

The response function decays as the power law t~^^'^. This makes a difference with 
the results obtained for the p-spherical model with p > 3, for which G{t,t') oc . It 
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implies that the memory ol the p = 2-model is too weak and the aging effects in the 
magnetisation are washed away quickly. 

The TRM decay for large times t, if the magnetic field has been applied during a finite 
interval [0,t^], behaves as m™^(t) oc r^/^ while if = A^t, > 0, rrJ^^{t) oc r^/^ 

One should note that many of the quantities calculated in this section depend on the 
initial conditions only via the term rji^jj^) = (s^(0))^. For these quantities the results would 
be unchanged if one took any random initial conditions such that rji^jj^) = 1, the overline 
indicating the average over the randomness in the starting configuration (conditional on 
knowing the matrix J). For example the distributions s^(0) = ±1 with equal probability 
and s^(0) ~ A^(0, 1) {i.e. zero mean Gaussian of unit variance) would have lead to the 
same results for such quantities. In addition the Gaussian initial conditions imply that 
Si(0) ~ A^(0, 1) by the rotational invariance of the Gaussian distribution. 

Finally, let us remark that expressions (n.5) and (IL14) for the correlation and response 
functions solve, without the use of any assumptions, the dynamical mean-field equation of 
the p-spherical model when p is set to 2. 

In the language of Refs. 18 and 19 the long time dynamics is described in terms of two 
functions, namely, the function X{C) that modifies the FDT 



and the 'triangular relation' / that relates any three correlation functions at three long times 



t > t" > t' . The fixed points of /, /(a, a) = a, separate different correlation scales. 

For the p = 2 model we know the exact solution for all times and we can then compare the 
long-time part of it with the results obtained in Ref. 18 with the help of these assumptions. 
In fact, at large times, the exact solution implies the existence of two correlation scales. 
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In the first one the correlation decays from I to qEA^ ^{C) = f and TTI and FDT are 
satisfied. This is the FDT scale. In the second one, the correlation decays from qEA to 0, 
X{C) ~ ^ and TTI and the FDT are violated. The result X = agrees with the 

predicted value X = [p — 2)(1 — q)/q of Ref. 18. 

As regards the triangular relation the situation is slightly more subtle. The two correla- 
tion scales are separated by the value C = qEA^ hence, the three fixed points of / are 0, qEA 
and 1. In general, just from its definition, between any two fixed points / must have the 
form 



with J a monotonous function to be determined by the model. There is going to be a function 
J for each correlation scale. 

In the FDT scale the correlation functions are homogeneous functions of time and hence 
one can always write a relation like (III. 24). One just uses the monotonicity of the correlation 
functions with respect to both times to invert the times in terms of the correlation functions: 



Ci2 = C(ti, t2) = C{h - h) ^h-h = C-i(Ci2), then Ci3 = C(C-i(Ci2) + C-i(C23)) and 
finally j(Ci2)=exp(C-HCi2)). 



In the second scale X = and the p = 2 long-time dynamic equations are identically 
satisfied without fixing the function j. The most one can say about the correlation functions 
with the formalism used in Refs. 18 and 19 is that an equation like Fq.(III.24) exists. 

It is interesting to note, however, that in the cases p > 3 the dynamical equations fix 
J to be the identity and then the triangular relation must be a product. They also imply 
C(t,t') = q hit') I hit) with hit) a monotonous function of time. In the case p = 2, where we 
know the exact solution, this simple scaling clearly does not hold. Instead one has something 
more complicated which, e.g. for T = 0, can be written as 



with h{t) = t and = 2^2 2/^/7(1 + uf'^- Only when A ~ one can write C(A) ~ 
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B. Staggered initial condition 



In this section we shall consider the case where T 7^ 0. When T = 0, if the system starts 
aligned in the direction of one of the eigenvalues of the interaction matrix then it is in a local 
energy minimum of the Hamiltonian and hence the dynamics is frozen. The computation of 
the functions F for the various staggered initial conditions is presented in Appendix A. 2. 

In the case a = 2, the function F leads to an exponential decay of the energy density 
towards £pq. When a ^ 2 the decay to £pq is as in the case of uniform initial conditions. 



The staggered magnetisation evolves in time as 



6[fi — a) 1 



;iii.26) 



Then, at large times 



a^2 



(3^^2(0) Stays zero Vt 

{S2{t)) y^qEAS2{0) 

(s^^a(t)) Stays zero Vt 

Thus, when the system starts from an initial configuration with a macroscopic condensation 
on the maximum eigenvalue, the averaged staggered magnetisation in this direction decays 
from its initial value to the asymptotic (macroscopic) one, with a weight ^J^ea- Instead, 
when the system starts from different staggered initial conditions, all the averaged staggered 
magnetisations are zero, asymptotically. 



The correlation with the initial conditions behaves as 
• If a = 2 



c(t,o)~y^ 



EA- 



;iii.27) 



• If a 7^ 2 



21 



C{t, 0) ~ ^e-(2-»)*y2^ (47r)4 (2t)t. (111.28) 

Hence, in contrast from the case of uniform initial conditions, the decay of C(t,0) is ex- 
ponential if the system starts condensed on the eigenvalues a 7^ 2, and C(t,0) tends to a 
constant if a = 2. 

In general for t ^ ^ 1 one finds using Eq.(.16) that 
• If a = 2 

+ tyj, tyj) ~ qEA- 

C(r + ,,, ,,) ^,,,2^ (yi±M) • . 2^2 [^) * . ,111.30) 

Hence for a ^ 2 the behaviour of the correlation function in this time regime exhibits 
explicitly the aging phenomenon. However, if a = 2 the system relaxes inside the equilibrium 
state of size qEA- In addition one can show that the full form of the correlation function 
(for sufficiently large t'), starting from the staggered initial condition a = 2, is precisely the 
same as the correlation function for uniform initial conditions restricted to the FDT regime 
(see Eq(III.12)). Hence in all large time regimes it satisfies time translation invariance. 

The two-replica overlap at = 0, Q{t,t), behaves as 
a = 2 Q{t,t) ^ qEA 

a ^2 Q{t, t) ~ const ■ q^^ ^ e'^^^'^^^ r^/^ (III.31) 

Hence, if the two replicas start from the 'equilibrium' initial condition, they evolve up to 
reaching a maximum distance of qEA^ the size of the equilibrium state. Instead if they start 
from any other staggered initial condition, they just separate completely, Q[t^t) decays 
exponentially with time. 
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(in.29) 



As far as the response function is concerned, it is easy to see that in the case where 
a ^ 2 the response behaves exactly as it does for the case of uniform initial conditions (for 
sufficiently large t'). For the case a = 2 one finds 



hence it clearly exhibits time translational invariance and exponential decay for small time 
differences. As was the story for the correlation function, it has for all large t', the same form 
as does the response function for uniform initial conditions restricted to the FDT regime. 
One may explicitly confirm that FDT is satisfied. 

Thus, any staggered initial condition with a ^ 2 fails to reach an equilibrium regime. 
The staggered initial condition a = 2 is the only one leading to equilibrium dynamics - 
exponential decays, time-translation invariance. 

Finally, let us mention that using the expression (.f3) derived in Appendix A for the 
auto-correlation function one can show that a delta-type divergence in the initial condition 
at /i = 2 is needed to ensure the equilibration of the system. 
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IV. TEMPERATURE VARIATION EXPERIMENTS 



In this Section we shall analyse the effects of temperature variations during the total 
waiting time in the asymptotic behaviour of the model. Since we are interested in non- 
equilibrium effects, namely aging effects, we shall let the system evolve from the uniform 
initial conditions of Section III, r][fi) = I, V/i. 

If the system starts at temperature T at t = and the temperature is changed to T + ST 
at a later time tj then the function F for this scenario is obtained as follows. For t < tj, it 
is given by the constant temperature result, FJ^^P(t) = Fx(t), as in Eq.(III.2). For t > tj, 
the computation of Appendix B gives 



In order to study this equation it is convenient to separate the 'asymptotic' factor exp(4t) 
(see Eq.(.7)) and to define the function 7: 







For times t such that t — t. <^ 1 




(IV.3) 



For times t such that t ^ tj ^ 1 Eq.(IV.l) can be approximated 



r 



2 ST 



df 7T+5T(t - f) 7T(t') 




(IV.4) 



For times t, ^ 1 and St = t — t. <^ 1: 
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^ £T{t) + {2 - T) 8T 8t . (IV.5) 

It is therefore clear from above that the energy density is continuous at the temperature jump 
(as one would expect on physical grounds). At first order in 8T the effect of the temperature 
change is linear in 8T . Fig. 3 shows the energy-density decay at constant temperatures 
T = 0.6 and T = 0.9. The dotted lines correspond to the equilibrium energy-densities 
£j = —0.7, —0.55, respectively. In addition, the curve corresponding to the temperature 
cycle r — )• r -|- 8T — > T with T = 0.6, 8T = 0.3, t^jj-^ = 50, ^^,2 = 100 is included. The 
effect of positive and negative temperature jumps can be seen in this curve. At t^jji the 
first perturbation is applied and the energy density grows to a value above the asymptotic 
energy-density £^_^^j, in an interval such that t — t^^ is short and then starts decaying 
towards the asymptotic value £^\gj'. A zoom in the figure would show that the perturbed 
curve for times bigger than t^^ is above the asymptotic energy £j\sj and below the curve 
associated to the relaxation at constant temperature T -\- 8T . It is an interesting feature 
of this system (and one that may be generic for non-equilibrium systems) that it reaches 
its final asymptotic energy density from above. Irrespective of the early time temperature 
pulse the time dependent energy density always becomes asymptotic to the energy density 
decay for a system that has been kept at constant (final) temperature. One can therefore say 
that a positive or negative temperature difference between the initial and final temperatures 
(irrespective of the duration of the pulse) has no effect on the ultimate rate at which the 
system reaches its final energy density. 

For very large times t as compared to the relaxation rate behaves as if the time-interval 
[0,tj] had not occurred: 

£j^mP(t) ~ £T+sT{t) = £T{t) + l^T. (IV.6) 
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Fig. 3. The thiner curves correspond to the energy-density decays at constant tem- 
perature T = 0.6,0.9. £j = —0.7,-0.55, respectively. The dotted hues describe the two 
equihbrium energies. The thicker hne corresponds to £{t) for the process T ^ T -\- ST — )• T, 
T = 0.6, 6T = 0.3, = 50, t^2 = 100. 



We now study the decay of the auto-correlation function when tj = t^. In this case, Eq. 
(II. 5) implies 



CJ™P(i.i„) = CT(i,i.) 



(IV.7) 



^ rJump(^) ' 

where Ct is the auto-correlation function at the 'initial' temperature T. The asymptotics for 
the auto-correlation function obtained in Section III therefore apply to the auto-correlation 
after a temperature jump, if one takes into account the presence of the multiplicative factor. 

• Short times after the temperature jump - FDT scale. 

At the beginning of the FDT scale t — t^, ^ 1, pJ^^P is given by Eq. (IV. 3) and 

cj^"^P(t, t^)^l-{T- 26T) [t - t^) + 0{{t - t^Y) . (IV.8) 



Instead, if t — t^, — )• oo, one can show that 
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(IV.9) 



where qEA{T) is the Edwards- Anderson order parameter associated to temperature T. 



• Large time-differences - Aging regime. 



In the limit of large times t ^ 



qEA{T + 6T) 

\J qEA{T) 



ST 
1 - T 



(IV.IO) 



It is clear from this expression that a positive (negative) temperature jump ST > 
(ST < 0) implies a smaller (bigger) correlation function and then a younger (older) 
system than the one associated to the initial temperature. The ST-eSect is linear at 
first order in ST. 

These results can be seen explicitly in Fig. 4, where we present a plot of the auto- 
correlation function vs. t [t = t — t^,), in a log-log scale, for tj = = 300 at constant 
temperature T = 0.3, 0.6, 0.9 and for the temperature jumps T = 0.6 and ST = 0, ±0.1, ±0.3. 
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Fig. 4. Auto-correlation decay in positive and negative temperature jump experiments. 
Log-log plot of C vs. T = t — t^jj. tj = t^jj = 300 in all the curves. The thiner curves 
are associated to constant temperatures T = 0.3, 0.6, 0.9 from top to bottom. The thicker 
curves correspond to T = 0.6 and 8T = —0.3, —0.1,0.1,0.3 from top to bottom, and hence 
to measuring temperatures 0.3,0.5,0.7,0.9 respectively. 

More interesting is to compare the decay of the correlation function after the temperature 
jump with that associated to the final temperature T -\- 8T , in the manner done in the 
experiments'*^. One sees that the curves after the jump get displaced in the direction of the 
older curves if 8T < or in the direction of younger curves if 8T > 0. This result is equivalent 
to the displacement of the maximum of the logarithmic derivative of the magnetisation decay 
observed experimentally'*. 

Nevertheless, this does not complete the understanding of the landscape of the model 
is concerned. The curves above are the analogous to the decay from the initial conditions 
C(t,0) elsewhere in the paper but where the initial conditions are those generated by the 
evolution of the system at constant temperature T for a time t^. 

If one performs the temperature jump at tj and then measures the C(t -|- 1^,, t^,) decay for 
t-w ^ tj then the rate of the correlation function decay can be said to depend on an effective 
age ta for the system at the final temperature. One would expect that ta ~ t^a -\- [t^ — tj) 
for some positive value a. In the case where a = then the effective age of the system is 
simply the time spent after the temperature change and hence the system has not benefitted 
from aging at the period [0,tj]. For t^, ^ one may verify that CJ*^*^P(tu, -|- r, t^,) ~ 
Cj+sjitw + T, tu,), i.e. the system has forgotten about the temperature jump. However if a 
is finite and ^ tj then C^^™^{tw + T^t^) ~ Cj^sTitjCi + {t^ — tj) + T,tja + {t^ — tj)), 
hence this asymptotic analysis is not sufficiently sensitive to reveal the effective age induced 
by waiting at a different temperature. Instead one is forced to consider the intermediate 
range where ~ 0{tj). This computation has been carried out numerically and the results 
are shown in Figs 5. a and 5.b. 
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One sees that for ST < (see Fig. 5. a.) the correlation function decay is very close to 
that of Cj+sjitw + T, tu,), that is a 1 and hence the effective age of the system is close to 
t^jj. In the case ST > one sees that the correlation decay is actually slower than that for 
the system at fixed temperature, implying that the effective age of the system (viewed at 
the final measurement temperature) is in fact greater than and hence a > I. 




0.1 1 10 100 1000 

t 



Fig. 5. a Negative temperature jump compared to constant temperature decays. C(t + 
tw) US T in a log-log scale. For the lower thin curve the temperature is constant, T = 0.6, 
and = 5. For the upper thin curve the temperature is constant, T = 0.6, and = 35, For 
the upper thick curve tj = 30 and = 35, T = 0.9, ST = —0.3, i.e. T = 0.9 — > T-\-ST = 0.6. 

It has been observed experimentally in TRM and out-of-phase susceptibility decays^ that 
a sample that has waited at a higher temperature [ST < 0) has an approximate effective age 
ta ^ t^jj — tj, that is the period spent at the higher initial temperature does not contribute 
to the effective age of the system. Fig. 5. a. shows the correlation decay for such a setting. 
It is apparent from the curves that this model does not capture this feature. The lack of 
chaoticity in the energy landscape of the p = 2 spherical spin glass means that it fails to 
reproduce this important feature of real spin glasses. 
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Fig. 5.b Positive temperature jump. C(t + t^^t^) us r in a log-log scale. For the thin 
continuous curve = 35 and T = 0.6. For the thick curve tj = 30 and = 35, T = 0.3, 
ST = 0.3, I.e. T = 0.3 ^ T + ST = 0.6. 

The response of the model to temperature variation is represented by the results for 
temperature jumps. The model we are here considering is not able to describe some of the, 
though controversial, most interesting features of spin-glass physics as is the asymmetry in 
the response to positive or negative temperature variations. 



One may also show that 

gjump^^ + t^, r + t^) = J'^l;^ Qt{t + t^, r + t^) (IV.ll) 

if the temperature jump is applied at the same time at which the noise is changed. Thus, 
for large times r ^ 1 

gjump^^ + t^, T + t^) ~ + + t^, r + t^) . (IV.12) 

V qEA{T) J 
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V. CONCLUSIONS 



In this paper we have examined the dynamics of the p = 2-spherical modeL In spite of 
having a very simple energy landscape, with only two equilibrium states and no signature 
of an exponential number of metastable states, its dynamical behaviour is far from being 
trivial. 

For almost any initial condition the system never reaches equilibrium, though its energy 
density asymptotically approaches that of equilibrium. Some of the experimental observa- 
tions are reproduced by the time-evolution of the auto-correlation function: a slow decay 
towards zero for widely separated times and aging effects represented by an explicit depen- 
dence on the waiting time. 

However, up to now there are no experimental results about the auto-correlation decay; 
the measurements are concentrated on the magnetisation decay. For this simple model the 
magnetisation depends on but it decays too fast, with an extra factor t~". It would we 
worth studying the effects of magnetic field variations in this solvable modeP^. 

A similar situation with aging effects in the correlation function and too fast a decay 
for the magnetisation has been obtained for 'unfrustrated' systems such as the XY modeP®. 
One may argue that the landscapes in which these models relax are far too simple and that 
this is the reason why the response function decays too fast. It is interesting to note that 
again, as it has already been noted in Ref. 28, the response appears to be more sensitive to 
the precise nature of the landscape than the auto-correlation function. 

The analysis allowed us to check, in every detail, the validity for this model of the weak- 
ergodicity breaking and weak long term memory (rather too weak in this case) hypotheses. 

The exact results demonstrate that, asymptotically, there are two time regimes. For 
close times compared to t', C decays from 1 to qEA- TTI is satisfied both by C and G and 
the FDT theorem holds. For well separated times C decays from qEA to 0, and TTI and the 
FDT theorem are violated. In the first regime C decays quickly while in the second regime 



31 



it decays slowly. 

As regards the function X[C] measuring the departure from the FDT theorem, one here 
finds X[C] = 1 if qEA 1^ C < 1 and X[C] = 0, if < C < qEA in good agreement with the 
extension of the results of Ref. 18 to the case p = 2. 

In the case of initial conditions leading the system to equilibrium dynamics [i.e. with 
initial macroscopic condensation on the eigenvector associated to the maximum eigenvalue 
/i = 2) we have been able to calculate the (time translationally invariant) autocorrelation 
and response functions. One finds that, somewhat surprisingly, these functions are identical 
to the corresponding functions in the FDT regime for initial conditions which ultimately 
lead to the aging phenomenon. It is not at all obvious that this should be the case and 
one is tempted to speculate on the possible generality of this phenomenon. It is possibly a 
consequence of a form of self similarity in the energy landscape; indeed we have also seen 
that all the initial conditions considered lead asymptotically to the same energy density^, 
adding further weight to this notion. 

We have also been able to examine the effects of temperature variations on the behaviour 
of the two-time autocorrelation function for this model. In real experimental situations one 
can assess a sort of effective age of a system that has experienced a temperature variation via 
TRM decay and a.c. susceptibility measurements ^ By analogy to the constant temperature 



^This is not expected to happen when p > 3. 

^To be more precise, in TRM experiments one can assign an effective age to a system that has 
undergone a negative cycling experiment. This is done by superposing the magnetisation decay 
curves after the cycle with a constant temperature decay for a smaller waiting time for all times 
explored experimentally; thus giving the effective waiting time. Nevertheless, this is not possible for 
positive cycling experiments, for which the decay at short times and at long times are considerably 
different, and the curves cannot be compared with constant temperature decays for any waiting 
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case, systems exhibiting a more rapid decay are said to be younger. A similar style analysis 
but using the correlation function to test the effective age has been carried out. Systems 
which have spent time at a higher temperature before switching to the final temperature do 
benefit from aging during time spent at higher temperature at variance with real spin glass 
experiments where the time spent at higher temperature does not contribute to the effective 
age of the system at the final temperature (if the temperature variation is big enough though 
all the temperatures are, of course, kept below the critical). 

Preliminary investigations would suggest a similar behaviour for the 0{N) model. How- 
ever, in this model one has a notion of space and we hope that a full study would allow us 
to develop an analytic picture of aging phenomena in terms of domain growth. Through 
this analysis we would like to attempt to make a connection between the droplet and the 
mean-field spin-glass models. 

In one respect the p = 2 spin glass shares a property of droplet models. The existence 
of two well defined ground states (at least at zero temperature), parallel to the eigenvector 
with maximal eigenvalue, implies that the system evolves via a competition of these two 
phases. Indeed, in common with this model, droplet models are unable to account for the 
experimental effects of temperature cycling unless the additional hypothesis of chaoticity is 
added by hand. 

It would be interesting to explicitly see how a complete 'aging' situation is settled when 
the parameter p, in the definition of the general p-spin spherical model, is raised from p = 2 
to a bigger value^^'^®. The solution we have here presented for the p = 2 model, is the exact 
solution of the full mean-field dynamical equations of Ref. 18 when p = 2. Hence, it can 
be used as the unperturbed solution in a perturbative analysis of those equations around 
p = 2. 

These final topics are currently under investigation^^. 



time. See Ref. 3 for a more detailed discussion on this point. 



33 



ACKNOWLEDGEMENTS 

We would like to thank G. Parisi for first drawing our attention to the interest of this 
model. During the course of this work we have benefited from useful discussions with J. -P. 
Bouchaud, J. Kurchan, M. Mezard and E. Vincent. L. E. C. wishes to thank the ITP at the 
University of California at Santa Barbara for hospitality. 



34 



Appendix A 



We here solve Eq.(II.7) at constant temperature T for F; we shall assume that we are in 
the region T < Tc {Tc = 1) - we shall a posteriori see that this guarantees the convergence 
of the expansions we use. The Laplace-transform of F is given by 

i-2r((^))' 

where r][fi) = (s^(to))^. Taking the inverse-transform 



T(s) = . (.11 

s — 2/i 



CO ^ k 

k=o i=l 

oo 



2{p - AO 



CO n n 1 /* 1 

+ Y.{2Tfk dX,p{X,)e''^' l[dX,p{X,)— dpp{p)^{p)— 

k=o 1=2 ^1^1 ~ ^1^1 ~ 

(•2) 

and re-summing the series 

V{t)= [dpp{p)rj{p)e'^' 



1 - Txip) 

+2T / .A piX) / dp pip) ,ip) ^ (.3) 



where 



X{p) ^PPjdX (.4) 

and PP indicates the Cauchy principal value. The initial condition Si(0) - or its 'staggered' 
distribution r][p) = (s^(0))^ - determines the time behaviour of F and, in particular, its 
asymptotic behaviour. 

If the interaction-matrix J belongs to the Gaussian ensemble, the density of its eigen- 
values is given by Eq.(IL8). For long times t the integrals over p are dominated by the 
maximum eigenvalue p = 2. Then, x(2) = 1 and 

rm~— I— (' de^er^{2-e)e-^'' + ^ T 6"^^* I dp ^^^^ ^^^U (.5) 

TT [ qEA Jo qjEA Jo J 2- e - p ) 
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A.l //-Uniform Initial Conditions 



If 7,(0) = 1, V/i, Eq.(.l)gi 



ives 



m = \;f.k^-^T\ (.6) 

fc=0 

for T < Tc = 1. Ik are the generalized Bessel functions. Its asymptotic behaviour is 

1 1 e'** 

^^^^ ~ 71^ (1 - r)2 (2t)^ ■ ^'^^ 



A. 2 Staggered initial condition 



If nip) = ^ then 



fc=o 



d\p{X) 



2{a - A) 



+A; / d\ p{\) e 



2Xt 



2(A - a) 



da p{a) 



and 



1 



T{t) = e2»* , ^— -T I dX '''''' -^^At 



2(A - a) 



1 



k-i ' 



:.8) 



(.9) 



l-rx(a) y a- A (l-rx(A))2- 

For long times the second term can be approximated. If Ji^ is in the Gaussian ensemble 
the saddle-point method implies 



V{t) ~ e 



2at 



1 T 



+ - e"* / de 



-2tt 



1 



l-Tx{a) TT Jo 2-a + e (l-rx(2-e))2 



(.10) 



• If a = 2 



Va = 2{t) 



.At 



(lEA 



(.11) 
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• If a 7^ 2 



o4t 



T 1 1 
71^ 4^ 2^ (2t)3/2 



A. 3 Computing the correlation function 

Once we have solved for T[t) we are left with the problem of computing the correlation 
function for non-zero waiting times t^. In order to do this we must calculate the term 



Wit,U) = / dr r(T) (( e^(*+2*»-2-) )) . (.13) 



Taking the Laplace transform with respect to yields 



W{t,s) = t{s){{-^)) (.14) 



and using Eq. (.1) yields 



W{t,s)=^-^^-^^ (.15) 

\ \ s — 2/i / / 

We now proceed as before, expanding the denominator then inverting the Laplace transform 
to obtain 



^ ^ i^-m-Txii^)) +i ^ ^ (A-/.)(A-/)(l-rx(/)) 

+/w ^^^.^^^^^^1^? ^L^^f , (.16) 



(^_^/)(i_rx(A))2 J ' 

this form is convenient for obtaining the relevant asymptotics in the case t ^ 
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Appendix B 



In this Appendix we compute the function T[t) for general cychc variations of the 
temperature. We start by obtaining the result for the single temperature jump case: 
T{t) = TOih - t) + (T + 8T)e{t - h). 

Let us define the function fit) 

m^{{e'''))= I diJ^ piij) e'^' (B.l) 



Starting from a uniform initial condition Eq.(II.7) implies, for the single temperature 
jump experiment, 

ri-jump(^^ = f{t) + 2T J df fit - t') r[o,,](t') + 2(r + 6T) J dt' fit - t') r[,,,o.](t') 

(B.2) 

where 

r[o,t,](t') = r(t)^(ti-t) (B.3) 

V^,^^^^it') = Vit)eit-tr) . (B.4) 

Solving the Laplace-transformed equation yields 

f = /» 2^r/(3)L[o,,](3) 

' 1 -2iT + ST)fis) 1 -2iT + ST)fis) ^ 

and taking the inverse-transform gives 

rJ^mP(t) = TT+5Tit) - 26T / dt' LT(t') TT+sHt - f) . (B.6) 

This solution should clearly yield LJ^^P(t) = Txit) for t < ti. This may be verified by 
noting that formally 

Lt = (i + 2r/)-V, (B.7) 
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where the products in the above expression are to be taken as the convolution of functions. 
Substituting this form for Tt+st into Eq. (B.6) yields the desired result. 

For completeness sake we mention that with a similar argument one can obtain the 
function F associated to n temperature jumps 

n 

T{t) = J2 T.-i e{t - u.,) e{u -t) + Tj{t - q . (b.s) 

8 = 1 

It is given by 

p.-jump ^ -2Y,{Tn - / dt' Trji - f) T^n-^^i^^V {^f) . (b.9) 

In particular, if n = 2^ To = T2 = T ^ Ti = T -\- ST and t > t2 we obtain the one cycle case: 

pi-cyc^^^ = Trit) + 26T /*' dt' Trit - f) v^-'^"^"^^ {t') (B.IO) 

■hi 

= VTit) + 2 ST r dt' Trit - t') VT+sT{t') 

rt2 rminiti ,t') 

-4 (STY / dt' Trit - t') / dt" VT+5T{t' - t") VT{t") . (B.ll) 

Jti Jo 
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